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Context. The observed spectral variation of HD 50138 has led different authors to classify it in a very wide range of spectral types 
and luminosity classes (from B5 to AO and III to la) and at different evolutionary stages as either HAeBe star or classical Be. 
Aims. Based on new high-resolution optical spectroscopic data from 1999 and 2007 associated to a photometric analysis, the aim of 
this work is to provide a deep spectroscopic description and a new set of parameters for this unclassified southern B[e] star and its 
interstellar extinction. 

Methods. From our high-resolution optical spectroscopic data separated by 8 years, we perform a detailed spectral description, 
presenting the variations seen and discussing their possible origin. We derive the interstellar extinction to HD 50138 by taking the 
influences of the circumstellar matter in the form of dust and an ionized disk into account. Based on photometric data from the 
literature and the new Hipparcos distance, we obtain a revised set of parameters for HD 50138. 

Results. Because of the spectral changes, we tentatively suggest that a new shell phase could have taken place prior to our observations 
in 2007. We find a color excess value of E(B-V) = 0.08 mag, and from the photometric analysis, we suggest that HD 50138 is a B6-7 
III-V star. A discussion of the different evolutionary scenarios is also provided. 
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1 . Introduction was cited by Merrill ( I1931I ) and Merrill & Burwell ( I1933l >. Since 

then, numerous papers have mainly considered this star as either 

Stars that present the B[e] phenomenon are known to form a a pre . main sequence star, more specifically as a Herbig Ae/Be 

■ heterogeneous group. This group is composed of objects at dif- star> or as a classical Be star . Jaschek et al. ( fl993T> have even 

C3 ; ferent evolutionary stages, such as high- and low-mass evolved conside red this star as a transition object between a classical Be 

stars, intermediate-mass pre-main sequence stars, and symbiotic and a B[e ] sta r. Later Lamers et al. ( fl998l and Zorec (TTWHl in- 

objects (Lamers et al.[1998). However, for more than 50% of the clude this star in the list of unc i assifled B [ e ] stars, 
confirmed galactic B[e] stars, the evolutionary stage is still un- 

known, so that they are gathered in the group of the unclassified The difficulty in obtaining the correct classification of 

B[e] stars. This problem is mainly caused by poor knowledge of HD 50138 18 mainl y caused b ^ stton^ spectral variability 

their physical parameters, especially their distances. ™ ln ^ ^^tf^^^^^.?f^^ 

In this paper, we present our study related to the southern 
B[e] star HD 50138 (V743 Mon, MWC158, IRAS 06491-0654). 



of different scales of spectral variability, from days to months, 

especially in the line profiles of Ha, Hei (A 5876), and Nai 

„, ,. . . , , ,. . , lines. Spectral variations have also been identified in the UV, 

The spectrum of this star was discussed tor the first time by , r , rTTT „ , TT „ | lr , - t „, 

vr-~ni -t „i /TTqtTI t a„ nrmmm . ~f b y the analysis of IUE spectra (Hutsemekers 1985 ). These spec- 



Merrill et al. (1925). Later, the presence of spectral variability , . 

' r r J tral variabilities have been explained by an outburst that possi- 

bly happened in 1978-1979 (Hutsemekers 1985) and by a shell 

* Based on observations (i) with the 1.52-m and 2.2-m telescopes phase in 1990-1991 (Andrillat & Houziaux [199TT ). On the other 

at the European Southern Observatory (La Silla, Chile), under agree- handi H D 50138 had presented small and non-periodic photo- 

ment with the Observatorio Nacional-MCT (Brazil) and (if) with the metrk variations that did not seem t0 be associated t0 the spec- 

Telescope Bernard Lyot, Observatory ol Pic du Midi (France). . , , , . , . , , , TT , ,| lnn a 

r J • tral changes related to this shell phase (Halbedel 1991 ). 

It is with great sadness that we have to report that, during the final 

stages of this paper, we had a deep loss when Francisco X. de Araujo In addition, studies based on polarimetry and spectro- 

passed away. polarimetry have identified an intrinsic polarization that seems 
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to be linked to nonspherical symmetry of matter around this ob- 
ject, probably a circumstellar disk (Vaidya et al. 1994, Bjorkman 
et al. [19981 Oudmaijer & Drew [l999l Harrington & Kuhn[2007]l. 
A disk-scattering effect associated to an outflow has also recently 
been suggested by Harrington & Kuhn (2009j>. 

On the other hand, Cidale et al. (2001) propose that this 
object could be a binary system. Based on spectro-astrometry, 
Baines et al. ([2006 ) suggest the same possibility, where the com- 
panion would be separated by 0.5" - 3.0". However, up to now, 
no direct evidence of binarity has been found. 

HD 50138 was further observed during the Hipparcos mis- 
sion, and its distance of 290+7 1 pc (Perryman et al. l 19971 > turned 
out with 500+1 50 pc to be almost doubled according to the new 
data reduction procedure performed by van Leeuwen (120071 1 . 
The newly determined distance and the large uncertainties in the 
previously published stellar classification attempt request and 
warrant detailed investigation and revision of the stellar param- 
eters of HD 50138. 

In this study, we present our new optical spectroscopic ob- 
servations and perform a photometric and spectroscopic analysis 
of HD 50138, aimed on the one hand at describing the observed 
spectral variations and, on the other, at better constraining the 
stellar parameters, needed for an improved discussion about the 
nature of this object. 

The paper has the following structure. In Sect.[2]we describe 
our observations. In Sect.[3] we present our results. In Sect. 13. II 
we describe our high-resolution spectra taken on different dates. 
In Sect. 13. 21 we derive from the analysis of our spectroscopic 
data and public photometric measurements, the insterstellar and 
circumstellar extinction, hence the stellar parameters and the 
spectral type of this object. In Sect.|4] we discuss the possible 
scenarios for explaining the nature of this curious star, and fi- 
nally in Sect.[5] we present our conclusions. 

2. Observations 

We obtained high-resolution optical spectra using the high- 
resolution Fiber-fed Extended Range Optical Spectrograph 
(FEROS) and the Narval spectro-polarimeter. The first FEROS 
spectrum was obtained on October 27, 1999, when the spectro- 
graph was attached to the 1 .52-m telescope, and the second was 
obtained on October 4, 2007, when it was attached to the 2.2- 
m telescope, both at the European Southern Observatory in La 
Silla (Chile). FEROS is a bench-mounted Echelle spectrograph 
with fibers, which covers a sky area of 2" of diameter, with a 
wavelength coverage from 3600 A to 9200 A and a spectral res- 
olution of R = 55 000 (in the region around 6000 A). We adopted 
its complete automatic online reduction, where the heliocentric 
correction is done. 

The spectrum of 1999 was obtained with an exposure time 
of 180 seconds and has S/N of approximately 80 in the 5500 A 
region. In 2007, we were able to take two consecutive spectra 
of the star, both with 180 seconds of exposure time. Since these 
spectra do not show significant differences, we added them up 
for a better S/N, which is around 250. 

Concerning the Narval data, they were obtained on March 
14, 2007. Narval is an Echelle spectro-polarimeter attached to 
the telescope Bernard Lyot at the observatory of Pic du Midi 
(France). For this study, we are only using the spectroscopic 
data that cover a sky area of 2.8" of diameter, with a wavelength 
range from 3750 A to 10500 A and a spectral resolution of R = 
80 000. We obtained 8 exposures of 300 seconds each. The S/N 
is around 360. We also adopted its complete automatic online re- 



duction; however, due to problems related to the merging of the 
spectral orders (especially in the region of the Balmer lines), we 
use this spectrum mainly for a qualitative comparison with the 
FEROS data, which are of better quality for our purposes. 

3. Results 

3. 1 . Spectral description of HD 501 38 

The high-resolution spectra of HD 50138 present lines from neu- 
tral and singly ionized elements. In addition to many emission 
lines with circumstellar origin, HD 50138 exhibits absorption 
lines, which are probably formed in the stellar photosphere. Our 
spectra were taken 8 years apart and strong spectral variations 
can be noted. Similar variations were previously reported in 
the literature and associated to shell phases and outburst events 
(Doa zan [19651 Huts emeke rs |T98"5l Andrillat & Houziaux [19911 
Bopp |T993l Pogodin [T9971 ). 

TableQ] lists all lines present in our FEROS spectra, along 
with their radial velocities (obtained considering the center of 
the lines and the laboratory wavelength of each transition), 
equivalent widths of emission and absorption components, and 
the possible identificatiorQ. For some lines there are more than 
one possible classification, however the radial velocity is derived 
assuming the first identification cited. There, "Uid" means that 
the line is unidentified. Because of uncertainties in the position 
of the underlying continuum, we estimate the errors of our mea- 
surements to be about 20% for the faint lines and about 10% for 
the strongest lines. Especially for the faint lines, there is also a 
significant difference between the radial velocities measured in 
1999 and in 2007. This difference probably comes not only from 
circumstellar changes, but also from the low S/N of our 1999 
data. The behavior of the lines from the main elements present 
in our spectra is described here. 

- Hydrogen 

In our FEROS spectra, we identify the presence of Balmer 
lines (see TableQ]). Four of them, HS, Hy, H/3, and Ha can be 
seen in Fig.Q] The profiles of these lines, except for Ha, are 
composed of a broad absorption component, probably of photo- 
spheric origin, and an emission and a narrow absorption com- 
ponents superimposed, which are probably formed in the cir- 
cumstellar medium. This kind of profile was already described 
by Houziaux (1960) from observations taken in the late fifties. 
The Ha line presents a double-peaked structure, where the V/R 
changes from 0.36 in 1999 to 0.88 in 2007. These values are 
different than those found by Halbedel (1991) and Pogodin 
(119971 ) of 0.6 and 0.5, respectively. However, the peak separa- 
tion presents small differences in both our spectra, ~ 160 km s _1 
in 1999 and 150 kms~' in 2007, in agreement with the value 
cited by Oudmaijer & Drew ( 1999} of ~ 160 km s -1 . 

Another significant variation is related to the radial veloci- 
ties of the narrow absorption components. In 1999, this value is 
almost constant for H6, Hy, and H/3, ~ 35 kms" 1 . However in 
2007, this value is quite different for each line, 95 kms -1 (H6), 
60 kms 1 (Hy) and 50 kms 1 (H/3). The Ha line's central ab- 
sorption presents a constant value of ~ 60 km s" 1 on both dates. 
This velocity is quite different from the velocity reported by 

1 In order to identify the lines, we used the line lists provided by 
Moore (1945), Thackeray < fl967t , and Landaberry et al. ( f200ll . We 
also looked up two sites on the web: NIST Atomic Spectra Database 
Lines Form (URL physics.nist.gov/cgi-bin/AtData/lines_form) and The 
Atomic Line List v2.04 (URL http://www.pa.uky.edu/~peter/atomic ). 
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Fig. 1. Balmer lines seen in the FEROS spectra of HD 50138. 
The solid line is related to the spectrum taken in 1999 and the 
dotted line to the 2007 one. 



Oudmaijer & Drew ( 1999 ); however, this can be related to differ- 
ent spectral resolutions. As cited by Pogodin i 19971 1, the Balmer 
lines also present symmetric wings, which show velocities that 
can reach ~ 1500 km s _1 , as can be seen in Fig.[T]and even better 
in the blow-up of the Ha wings in Fig. [2] However, these high 
velocities seem to come from electron scattering, since no other 
element identified in our data present such characteristic. 

We identify Paschen 22 to Paschen 12 lines in our spectra. 
We show one of these lines, Paschen 19, in Fig. [3] They dis- 
play strong variations in our FEROS and Narval data. In the 
1999 FEROS data, we detect a single-peaked profile with a blue 
shoulder. This shoulder becomes more prominent, probably be- 
cause of the better S/N, in the Narval data. However, with the 
2007 FEROS data, we clearly have a double-peaked profile with 
a peak separation of 190 kms" 1 . Assuming the shoulder previ- 
ously seen in 1999 and in the Narval data as another peak, we 
have a variation of V/R in the 2007 FEROS data. A similar vari- 
ation was reported by Andrillat & Houziaux ( 119911 ), which was 
associated to an outburst. The wings of Paschen lines are ex- 
tended in all spectra to ~ 300 km s _1 . However, as seen in Tabl^U 
we believe the radial velocities obtained for the Paschen lines are 
affected by some problems with the precision of the wavelength 
calibration of FEROS, beyond 8300 A. 

- Iron 

Iron, as can be seen in TableQ] is by far the element with the 
most of lines identified in our spectra, with all of them from Fen. 
We have identified permitted and forbidden lines. The latter are 
represented by narrow single-peaked emission lines from multi- 
plets 3, 4, 6, 7, 17, 18, 19, 20, 21, and 34. These forbidden lines 
do not present any noticeable difference in all spectra. 

We have identified permitted lines from the multiplets 20, 2 1 , 
25, 26, 27, 28, 30, 37, 38, 41, 42, 43, 49, 55, 73, 74, 167, 190, 




-1000 
Heliocentric Velocity (km/s) 



Fig. 2. Ha wings seen in the FEROS spectra of HD 50138. The 
different line styles have the same meaning as in Fig. 1. The 
normalized continuum is also presented for better visualization. 
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Fig. 3. Variations in the line profile of Paschen 19 at 8413A. 



200, and 203. Unlike Jaschek & Andrillat ( 1998), where the Fen 
lines show absorption (mainly) or emission profiles, our spec- 
trum presents a special kind of profile for each date (see Fig.|4|). 
In the 1999 FEROS spectrum, these lines have a P-Cygni profile, 
where the absorption component is made up of two components, 
with the bluer one being more intense. In the 2007 FEROS spec- 
trum, these lines have a different profile. There is an emission 
profile associated to a strong central absorption that also has two 
components, where the bluer one is also more intense. These ab- 
sorption components are shifted to higher velocities compared 
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Fig. 4. Variations in the line profiles of the Fen multiplet42. The Fig. 5. Profiles of the permitted lines of Oi and their variation, 
bottom line corresponds to the line at 4923A, the middle one at The bottom line corresponds to the triplet at 7774A and the top 



50 18 A, and the top one at 5 169 A. 



one at 8446A. 



to 1999. On the other hand, the Narval data show a sort of tran- 
sition between the P-Cygni profile and the profile seen in the 
2007 FEROS spectrum. The absorption components are not as 
intense as seen in the other spectra; however, three of them can 
be seen, where the red one, contrary to the FEROS data, is the 
most intense. 

- Oxygen 

Permitted and forbidden lines of neutral oxygen are seen in 
the HD 50138 spectra. No On lines were identified, indicating 
that the emission region is not highly ionized or that the higher 
ionized regions are rather small in volume, supressing the effi- 
cient production of observable On emission and absorption. The 
forbidden lines present narrow emission profiles and do not show 
any significant differences in all the spectra, with an almost con- 
stant radial velocity around 40 km s _1 (Table[TJ. All three lines 
from the multiplet 1 are seen, different than reported by Andrillat 
& Houziaux (1972). In addition, a contamination by airglows, 
based on the comparison with spectra from the sky taken at the 
same time, is discarded. 

The permitted lines present a complex profile and show a 
strong variation in the different spectra (see Fig. [5}. The line 
centered at 7774A is in fact a triplet, thus the profile seen in 
our data has a combination of possible emission and absorp- 
tion components from each line. However, the profile of this 
triplet changes strongly at different dates, as already mentioned 
by Jaschek & Andrillat (fl998l . Unlike the 1999 data, the 2007 
FEROS data show the absorption components clearly superim- 
posed over a strong emission component, similar to the Fen 
lines. In the Narval spectrum, this line has a transitional profile 
between those detected in the different FEROS spectra, espe- 
cially in the emission component. 

Another oxygen line present in our spectra is located at 
8446A. This line presents a strong emission in the FEROS 
spectrum of 1999 and in the Narval data; however, this profile 



changes completely in the FEROS spectrum of 2007, because it 
is less intense and has multiple peaks. 

- Helium 

HD 50138 presents several Hei lines from multiplets 2, 10, 
11, 12, 14, 16, 18, 20, 22, 46, 47, 51, 53, and 55. Line pro- 
file variations are clearly seen in the different spectra (Fig.|6]l. 
In 1999, we can see an inverse P-Cygni profile for some lines, 
such as 5876A. The separation seen between the emission and 
the absorption components agrees with the value cited by Bopp 
( 119931 ), ~ 200 km s _1 . However, in 2007, profiles similar as seen 
for the Fen lines are formed, where the absorption component 
is narrower than in 1999. Both profiles present wings with sim- 
ilar extension velocities, as can be seen in Fig. [6] Similar varia- 
tions were described by Bopp dl993l l, who reported an inverse 
P-Cygni profile, a few months after the outburst in 1991, and 
also by Pogodin dl9971 l, who cited that the profile of the line at 
5876A, changes from night to night. 

Our spectra do not show any Hen line. The tentative de- 
tection of the Hen line at 10123 A, as reported by Jaschek & 
Andrillat ( 1998 ), could not be confirmed by our Narval data. We 
can therefore claim that this object is not very hot. 

- Magnesium 

Only one line from Mgi, at 63 18 A, is identified in our spec- 
tra. On the other hand, few lines from Mgn are seen in our data, 
with the most prominent one, at 4481 A. This line mainly has 
photospheric origin, and as can be seen in Fig.|7j it also presents 
a noticeable variation in our FEROS data taken with 8 years of 
difference. 

- Silicon 

In the FEROS spectrum of 1999, Sin lines have pure absorp- 
tion profiles. On the other hand, in March 2007, the lines of the 
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Fig. 6. Line profile variation of the Hei line at 5876A. The solid Fig. 8. The variation of the line profile of Sin at 6347A. 
line is related to the spectrum taken in 1999 and the dotted line 
to the 2007 one. 
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Fig. 7. As in Fig.|6]but for the Mgn line at 448 1 A. 



multiplet 2 (AA 6347,6371) developed into an inverse P-Cygni 
profile, which is even more pronounced in October 2007, as can 
be seen in Fig.[8]for one of these lines. However, the lines of the 
multiplet 3 (AA 4128,4131) still show pure absorption profiles at 
these dates (see Fig. [9} and might thus be considered as having a 
pure photospheric origin. 



Fig.9. As in Fig.|6]but for the Sin lines at4128A and4131A. 



interstellar contributions. However, it is impossible for us to dis- 
entangle each component, which makes it impossible to derive a 
distance to HD 50138 from these lines. These lines also present 
variability from 1999 to 2007 (Fig.fTOli. because it is almost a 
P-Cygni profile, with three absorption components in 1999 and 
a broad emission component superimposed by three absorption 
components in 2007. 



- Sodium 



- Other lines 



Lines of Nai are present in our spectra at 5890A and 5895A. Permitted lines from Tin, Can, Cm, Cn, Sn, and Ni are also 
These lines are formed by a combination of circumstellar and identified in our spectra (see Table[TJ. These lines also show sig- 
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Fig. 10. Different line profiles of the Nai lines at 5890A and 
5896A. The normalized continuum is also presented for each 
date. 



nificant differences from 1999 to 2007, following the same be- 
havior as described for the Fen lines. 



3.2. Determination of the extinction and spectral 
classification 

The determination of the stellar parameters of HD 50138 is diffi- 
cult, since both its spectral lines and its continuum show strong 
variations. It is, therefore, not surprising that the range of spec- 
tral types (from earlier than B5 down to AO) and luminosity 
classes (I-V) found in the literature is rather wide, because we 
can expect that these determinations depend on the method (or 
lines) used for the classification, as well as on the period at which 
the observations were obtained. 

To constrain the spectral classification of HD 50138, we 
made use of both our own sets of high-resolution optical spectra 
and public photometric data. We collected all available UBVRI 
band photometric data. These observations cover a time interval 
of almost 25 years and are displayed in Fig.Qj] The references 
for the data are listed in Table[2] 

Unfortunately, hardly any photometric data could be found 
before and around the time of the outburst, which was reported 
to have happened in 1978/79 (Hutsemekers |1985l l. The conse- 
quences of this outburst are, however, clearly visible in the light 
curves displayed inFig.QTl afading of the star in the £/BVbands 
by about 0.2 mag, and a recovery to the old brightness after about 
5 years. This trend in photometric data might be explained by 
some large mass ejection during the outburst, hiding the stel- 
lar continuum inside a dense shell or ring of (partly) optically 
thick material, which then turns transparent during expansion, 
allowing the underlying star to be seen again. On the other hand, 
the shell phase seen in the optical spectra in 1990/91 (Andrillat 
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Fig. 11. Time variation of photometric observations in the 
UBVRI bands collected from the literature (see Table|2]i. Also 
shown are the supposed positions of the outburst and the shell 
phase reported in the literature. The dashed lines in the B and V 
bands indicate the mean values before and after the outburst. The 
dotted bar in the B and V light-curves refer to the observations 
of Alvarez & Schuster GMT). 



Table 2. References to the UBVRI observations of HD 50138 
shown in Fig.fTTI 



Julian Dates 


Reference 


[2400000+] 





39916-40464 
40556-41224 
43350 - 44200 
44589 - 44596 
44934 - 45809 
46392-48331 
48699 - 48707 



Penston ( fT973l 
Alvarez & Schu ster <fT98ll 
de Winter et al. d2001b 
Kilkenny et al. ( 1985) 

Halbedel il99l\ 

de Winter et al. ( t200TT > 



& Houziaux 1991) remained invisible in the photometric data, 
as cited by Halbedel (1991). This might indicate that the shell 
phase was caused by mass ejection with much less mass loss. 

Interestingly, after the recovery of the star, it seems to be 
even brighter than it was before the outburst, as can be seen by 
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the comparison of the mean values before and after the outburst 
(Fig.fTTTi. We also notice that the mean value before is based 
only on two datapoints, but even considering the large errorbars 
of these observations, it seems that the star appears brighter af- 
ter ~ 1984. In addition, Alvarez & Schuster ( 119811 ) observed 
HD 50138 in a period close to the outburst, but without provid- 
ing either the exact dates or the observed magnitudes. However, 
these authors claimed that during their observations, HD 50138 
was not variable, compared to previous observations, meaning 
that AB and AV < 0.02 mag, according to their definition of 
variability. 

This apparent brightening after ~ 1984 might be caused by 
some contribution of the circumstellar material ejected during 
the outburst and after adding flux at all bands. Even though 
HD 50138 is known to possess a strong infrared excess emis- 
sion due to circumstellar dust (Allen [T973l , this dust will hardly 
contribute at optical wavelengths, i.e., in the UBV bands. On 
the other hand, based on spectropolarimetric observations per- 
formed in 1995 by Bjorkman et al. ( fT998j l, HD 50138 seems to 
have an almost edge-on ionized Be-like disk. Such a disk, which 
might have been formed from the outburst material in a similar 
way to how disks are formed around classical Be stars, can well 
add flux to the optical continuum in the form of free-free and 
especially free-bound emission. 

Whether the additional continuum emission from the ionized 
circumstellar material/disk will hamper proper spectral classifi- 
cation of HD 50138, or whether we even have to account for cir- 
cumstellar dust along the line of sight towards the star, acting as 
a further (i.e. circumstellar) extinction source, will be discussed 
and investigated in detail in Sects. [3~l2~3"l and [3~l2~2l respectively. 

3.2.1. The interstellar extinction 

First, we estimate the total extinction towards HD 50138, i.e., 
neglect any circumstellar contribution and thus derive an upper 
limit on the interstellar one. For this, we make use of the two 
color indices (U - B) and (B - V). Their variations with time are 
plotted in Fig. [131 

Obviously, (U - B) is quite variable, while (B - V) remains 
rather stable, even during the outburst phase. Despite the vari- 
ability in ( U - B), we can obtain a reasonable range for the possi- 
ble stellar classification from these observed color indices, mak- 
ing use of the following relations for the interstellar extinction 




-0.6 



-0.15 



(B-V) 



-0.1 



Fig. 13. Constraining the range in spectral type for different lu- 
minosity classes. The solid lines correspond to the upper and 
lower limits derived from observations. 



determinations, 



E(B -V) = (B-V)-(B- V) 
E(B - V) 



(1) 



J^E(U-B) = J-[(U-B)-(U-B)o} . (2) 



The second relation follows from Leitherer & Wolf (1984), and 
the parameters (U - B)q and (B - V)o refer to the intrinsic color 
indices. By combining these two equations, we obtain for each 
observed set of (U - B) and (B - V) colors a reddening indepen- 
dent equation of the form 



(C/-B) = 0.77{(B- V)o-*}, 
relating the two intrinsic colors, with 



x = (B-V) 



(U-B) 
0.77 



(3) 



(4) 



We calculate the parameter x, hence the relation Eq. ©, for 
each observation. Then we search for the two boundary cases, 
i.e., the highest and lowest values for x, and plot the resulting 
range in relation Eq. (01 in Fig. [13] Next, we look up the tables 
of Schmidt-Kaler (1982) for the intrinsic colors of stars in the 
spectral range B3 to Al. The luminosity class of HD 50138 is 
not well known, but previous classifications found in the litera- 
ture tend towards luminosity classes III-V (e.g., Houziaux 1960; 
Houziaux & Andrillat |1976l Fremat et al. 2006). Based on the 
Hipparcos distance, we can already exclude a supergiant classi- 
fication for HD 50138, since it would deliver an unrealistically 
high extinction value. We thus restrict it to luminosity classes II 
to V. 

The behavior of the intrinsic colors from Schmidt-Kaler 
(119821 ) is included in Fig. [13] From the overlap with the ob- 
servational ranges we find some preliminary sets, as listed in 
Table[3] of MK types, effective temperatures and color excess 
values, which are upper limits to the real interstellar extinction. 
The assignment of the temperature ranges is made based on the 
tables of Flower ([T996j ). 
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Table 3. Possible stellar classifications of HD 50138 under the 
assumption of pure interstellar extinction along the line of sight. 



L.C. 


Sp.Type 


T eff [K] 


E(B - V) mllx [mag\ 


V 


B5.5-7.5 


13 300 ± 900 


0.15+0.03 


III/IV 


B5.5-8.0 


13 100 ± 1 100 


0.15 ±0.03 


II 


B7.5-8.5 


1 1 600 ± 400 


0.11 +0.02 



A rather low value for the interstellar extinction has al- 
ready been suggested by Houziaux & Andrillat (1976) and 
Hutsemekers (1985 ) based on the weakness of the interstellar 
feature at /12175A and on the absence of interstellar lines in 
the IUE spectrum. To split the total extinction value into its 
interstellar and (possible) circumstellar contributions we, there- 
fore, search for further extinction indicators in our spectra. One 
such indicator makes use of diffuse interstellar bands (DIBs). 
In our high-resolution spectra, we found one reasonably good 
DIB at /15780 A, from which the extinction can be derived (see 
Herbig [T993l l. This feature is present in both the 1999 and the 
2007 FEROS spectra, but because of the poor S/N of the 1999 
data we can rely only on our 2007 data. From the equivalent 
width of 0.04 A of this DIB, we derive an extinction value of 
E(B - V) = 0.08 + 0.01 mag. 

As a further check we investigate the interstellar extinction 
within the galactic plane in the region around our target. It turns 
out that HD 50138 is located within a region of rather low ex- 
tinction (see Neckel & Klare |1980l Chen et al. |19891 Arenou et 
al. [19921 . Using the model developed by Hakkila et al. (1997) 
for the large-scale visual interstellar extinction, assuming Ry = 
A V /E(B-V) = 3.1, we have obtained Ay = 0.22 + 0.22 mag, de- 
livering E(B - V) = 0.07 + 0.07 mag. This value agrees with the 
value we found from the DIB. However, because of higher un- 
certainty (and definitely larger error) with this method, the value 
found from the DIB of E(B - V) = 0.08 + 0.01 mag seems to be 
reasonably accurate, leading us to use this value later. Compared 
to the (upper limit) values we found from the color indices, the 
contribution of the interstellar extinction towards HD 50138 is 
thus only half of it, indicating that the other part must be cir- 
cumstellar in nature. 
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Fig. 14. Intrinsic colors as a function of dust optical depth. 
Individual curves represent calculations performed for a mixture 
of silicates and amophous carbon grains for either single (small- 
est and largest, as cited in the text) grain sizes or an MRN grain 
size distribution. The dotted lines indicate theoretically expected 
intrinsic colors for different spectral types. 



The dust extinction parameter, Af??L, can be calculated from 



3.2.2. The influence of circumstellar dust 

The presence of a strong infrared excess implies that HD 50138 
must be surrounded by circumstellar dust (Allen 1973 ). How this 
dust is distributed, i.e., whether it is situated within a disk or 
shell, within or outside the line of sight, is not known. 

The circumstellar dust, like the interstellar one, absorbs and 
scatters the light from the star, so that we can define a circum- 
stellar extinction parameter, A^JL The intrinsic magnitude of 
the star in each photometric band, UBVo, is then obtained from 
the observed one, corrected for the amounts of interstellar and 
circumstellar extinction, i.e, 



UBVo = UBV ohs 



1UBV 



I ISM 



v /ISM 



X UBV ■ 



(5) 



with the parameters (Aubv/Av)ism resulting from the mean in- 
terstellar extinction curve (see e.g. Cardelli et al. 119891 Mathis 
1990). A star that suffers from additional circumstellar extinc- 
tion thus appears fainter. 



* dust 



1 



0.41n(10) of 



UBV dust 



(6) 



with <x ext as the extinction coefficient of the dust within the corre- 
sponding band and T^, ust as the optical depth of the dust at visual 
wavelengths. The extinction coefficients depend on the size of 
the grains and their chemical composition. Since the dust com- 
position along the line of sight is not known, we allow for dif- 
ferent sizes of the (spherically) test grains, which we consider as 
composed either of amorphous carbon or silicates. The range of 
sizes used extends from 7.5 x 10~ 3 fim up to ~ 1 /im. In addition 
to single grains of different sizes, we consider the well-known 
MRN grain size distribution (Mathis et al. 119771 . 

For each grain size and species, the dust extinction depends 
on only one free parameter, which is T?, ust . Expressing the set of 
Eqs. (0 by the intrinsic color indices (U - B)o and (B - V)o, we 
end up with two equations 



(U-B) = /(r d v ust ), 



(B - V) = f(4 ust ) 



(7) 
(8) 
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These calculated intrinsic colors for different dust optical depths 
can then be compared to theoretically expected color indices 
of stars with different spectral type and luminosity class (e.g. 
Schmidt-Kaler 1982 ), in order to constrain the range of possible 
spectral classifications. 

The results are shown in Fig. [14] for stars with luminosity 
class V, but the results for the other luminosity classes are very 
similar. Plotted are the dependences of the two intrinsic color 
indices on r^ ust computed for the smallest and largest grains, as 
well as for the MRN grain size distribution, using a mixture of 
silicates and amorphous carbon. The theoretically expected val- 
ues for the different spectral types are also shown. Obviously, an 
agreement between the computed and theoretical color indices 
(delivering the same dust optical depth) is achieved only for stars 
with spectral type between B6 and B7. Very similar results are 
found for the luminosity classes III/IV, while for stars with lu- 
minosity class II we find a possible spectral range of B8-8.5. 

Interestingly, pure large grains are not able to account for the 
circumstellar extinction. Instead, the circumstellar dust along the 
line of sight must consist of predominantly small grains, which 
can be located in an optically thin dust sphere or shell, in agree- 
ment with the results of Bjorkman et al. dl998l l. However, other 
scenarios, like a dusty disk where the dust is not in our line of 
sight or is optically thin, since it is formed far from the star, can- 
not be discarded. An interferometric analysis will probably an- 
swer this question about the circumstellar dust geometry (Borges 
Fernandes et al., in preparation). 

3.2.3. The influence of the ionized circumstellar disk 

The contribution of an ionized envelope or disk to the total con- 
tinuum emission of early-type stars was the subject of many de- 
tailed investigations during the past years. For the case of classi- 
cal Be stars, for instance, Zorec & Briot ( 119911 ) found a relation 
between the excess emission in the V band (defined as the dif- 
ference between the V band magnitudes measured during the Be 
phase and the non-Be phase) and the effective temperature of the 
star. Their relation indicates a stronger circumstellar contribution 
for hotter (i.e., earlier) Be stars, while late-type Be stars hardly 
show any contamination of their V band fluxes by their circum- 
stellar ionized material. A similar trend was found by detailed 
numerical studies of Stee & Bittar ( 120011 ), who found circum- 
stellar contributions to the continuum in the B and V bands for a 
B5e star of 0.57% and 6.01%, respectively, while the contribu- 
tions are found to be much lower for later types. But not only are 
classical Be stars influenced by their ionized circumstellar ma- 
terial, B-type supergiants can also have a significant wind con- 
tribution at all optical bands, as recently shown by Kraus et al. 
(120081120091 . 

B[e] stars are generally believed to have disks with higher 
densities than classical Be stars. In this case, an influence of the 
free-free and at optical wavelengths, especially of the free-bound 
continuum emission, might thus be expected. But in contrast to 
circumstellar dust, which acts at optical wavelengths as a pure 
absorber that dims the stellar light, the free-free and free-bound 
processes act not only as an absorber, but at the same time as 
an additional emission component, which adds more flux at a 
given wavelength than it absorbs from the underlying stellar flux. 
A star with an ionized disk thus usually appears brighter at all 
wavelengths compared to a star with no ionized disk (see Kraus 
etal. l20U8l 120091 . 

We check the influence of such an ionized circumstellar disk 
on the proper spectral type determination by calculating the 
emission of free-free and free-bound processes in circumstellar 



disks of different densities. The shape of the continuum emis- 
sion of the ionized gas does not severely depend on the geome- 
try of the ionized material. For instance, a spherically symmetric 
wind and a geometrically flat disk result in the same shape of 
the wavelength dependent free-free and free-bound continuum 
emission. This indicates that the identical continuum flux dis- 
tribution can be obtained from quite different geometrical sce- 
narios. In addition, the free-free and free-bound continuum at 
optical wavelengths is generated in the close vicinity of the star, 
typically within 2 - 5 R, (Kraus et al.[2008 ). It can thus not offer 
any insight into the global disk density distribution and geome- 
try. 

Since our aim is only to obtain proper stellar parameters, 
but not a global description of the detailed geometry and den- 
sity variations within the highly dynamical and non-spherically 
symmetric circumstellar material, for which we definitely have 
too little observational information at hand, we study the influ- 
ence of the ionized disk emission on the color indices in a very 
simplified way. 

From the literature values (see Fig.fTTI) we computed the 
pre-outburst average observed magnitudes in the UBV bands, 
deredden them with the interstellar extinction value derived in 
Sect. l3.2.Tl and convert them into fluxes. These fluxes corre- 
spond to the star plus disk system. Then, individual disk fluxes 
in the UBV bands were calculated. For this, we adopted a very 
simple disk model whose surface density distribution can be de- 
scribed by 



S(r) = 1.0* — . 

r 



(9) 



Such a surface density distribution follows, e.g., for the case of 
an outflowing disk forming wind, which might be an appropriate 
model. In that case, the surface density at the inner edge of the 
disk is given by (see Borges Fernandes et al. 120071) 



Ma 



(10) 



where M describes the (here constant) mass loss rate over the 
disk-forming wind region, a is the disk opening angle, and u the 
disk outflow velocity. 

The disk fluxes (in the form of free-free and free-bound 
emission) in the UBV bands, resulting from disks with increas- 
ing values of Eo, are then subtracted from the observed and inter- 
stellar extinction corrected ones. The resulting pure stellar fluxes 
are converted back into magnitudes and the color indices are de- 
rived. The results are shown in Fig. [15] Obviously, for low sur- 
face densities, the contribution of the disk is negligible, while its 
importance drastically increases with increasing surface density 
(top panel of Fig. [TBI. 

Because of the increasing disk flux in all three bands with 
respect to the stellar flux, the star must be intrinsically fainter 
in all three bands (see Fig. [TBI, especially influencing the (B-V) 
color index (lower panel) more than the (U-B) (middle panel). 
From a comparison of the derived color indices to the expected 
intrinsic colors of Schmidt-Kaler (1982) for stars of luminosity 
class V and different spectral types, we find a rather narrow range 
of valid stellar classifications, as indicated by the box extending 
over the complete plot in Fig.[TBJ While for luminosity classes 
III-V the range in spectral type is B6-7, we find a spectral type 
of B8 for luminosity class II. These ranges are in fairly good 
agreement with those found previously. This is not surprising, 
because the amount of circumstellar extinction is rather small, 
so that no big changes in intrinsic colors can be expected. 
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Fig. 15. Increase in disk flux with respect to the stellar flux (top 
panel) and intrinsic colors (mid and bottom panel) with disk sur- 
face density. The dotted lines indicate theoretically expected in- 
trinsic colors for different spectral types and a luminosity class 
V. The box extending over all three panels gives the range in 
possible classifications and the corresponding needed disk flux 
in all three bands. 



We would like to emphasize that the surface density chosen 
for the x-axis to plot our results in Fig. Q3] is not the most rel- 
evant parameter. Instead, the influence of the disk as shown in 
the top panel is the crucial parameter here. Such fractions of 
the free-free and free-bound continuum can also be achieved 
with quite different geometrical models and density distribu- 
tions. Nevertheless, the range in values for "Eo found from our 
analysis can easily be achieved with typical values for B[e] stars, 
e.g., disk opening angles in the range of 5 - 30 degrees, disk out- 
flow velocities of 10 - 100 kms _I , and mass loss rates typically 
of 10~ 8 to 10~ 5 M yr . However, the exact value for each pa- 
rameter cannot be provided by our research due to the lack of 
reliable information concerning the geometry and dynamics of 
such a possible pre-outburst ionized disk, and the modeling pro- 
vided here cannot be regarded as the only valid answer. It was 
only used to qualitatively discuss the influence of the ionized 



material to the color indices of the star in order to constrain the 
spectral classification of HD 50138. How the ionized material is 
really distributed around the star cannot be derived from a pure 
free-free and free-bound emission calculation. 

From our analysis of the circumstellar extinction contribu- 
tions of either the dust or the ionized disk, we find possible clas- 
sifications of HD 50138 as either a B6-7III-V star with T eS = 
13200 + 500 K, or a B8-8.5 II star with r eff = 11300 + 300 K. As 
a further check of these derived classifications, we derived the 
ratio of the equivalent widths of the most plausible photospheric 
Sin lines in our spectra. The ratios of the 4131 A/4128 A and the 
5056 A/5041 A equivalent widths are sensitive to temperature 
and surface gravity and have been computed from the line identi- 
fication tables obtained with the code S YNSPEC (see Hubeny & 
Lanz 2000) based on Kurucz ( 119791 ) model atmospheres in local 
thermodynamical equilibrium and Kurucz (1993) line lists. At 
temperatures between 10000K and 12 000K and surface grav- 
ities according to luminosity class II, these ratios show a mild 
increase with temperature from 1.7 to 1.8 (4131 A/4128 A) and 
from 2.1 to 2.3 (5056 A/5041 A), while the observed values, 
based on TableQ] are 1.41 + 0.19 and 1.82 + 0.26 in 1999 and 
1.08 + 0.15 and 1.50 + 0.21 in 2007. From this discrepancy be- 
tween observed and calculated equivalent width ratios, a B8-8.5 
II scenario for HD 50138 can be excluded. For higher temper- 
atures and especially higher surface gravities, the variations in 
the equivalent width ratios are no longer strictly correlated, but 
an overlap in the theoretical and observational values for stars 
of spectral types in the range B6-7 III-V occurs, so that these 
equivalent width ratios can be used at least as a qualitative check 
for consistency. From the photometric results and the qualita- 
tive agreement with the equivalent width ratios of Sin lines, we 
can claim that the final spectral classification of HD 50138 is a 
B6-7 III-V star. It is important to cite that our results agree with 
the classification obtained by Zorec et al. (1998) and Cidale et 
al.j2001 ) using the BCD (Barbier-Chalonge-Divan) spectropho- 
tometric system, based on the study of the Balmer discontinuity. 

To finish our stellar parameter study, we want to deter- 
mine the possible ranges in stellar luminosities. Our ionized 
disk model implies a pure stellar V band flux from 6.49 mag 
to 6.71 mag for the lowest and highest possible disk contribu- 
tions of 7% and 20%, respectively, in the range of possible MK 
types. We thus conclude that the star has an intrinsic V band 
flux of 6.60 + 0.11 mag. This result has been obtained from con- 
sidering the pre-outburst photometry only. We also tried it with 
the post-outburst data. Thanks to the high variability of the data, 
the results still agree with a B6-7 star but have a higher uncer- 
tainty. Using the bolometric corrections for B6-7 III-V stars from 
Flower (1996) and the Hipparcos distances towards HD 50138 
of 500 ± 150pc, we can finally calculate the stellar luminosity to 
log(L«/L ) = 3.06 + 0.27. 

We did not attempt to derive the luminosities for the case of 
the circumstellar dust extinction, because (i) based on the ob- 
served ionized disk around HD 50138, this possibility seems to 
be the more realistic one, and (ii) we currently have no clear 
indication of the possible dust composition and grain size distri- 
bution at hand. Nevertheless, assuming the range in dust optical 
dephts obtained by our analysis for the MRN grain size distri- 
bution (see Fig.[T4b. the range of stellar luminosities corrected 
by the circumstellar dust extinction is in a fairly good agreement 
with the above derived luminosity range. 
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4. Discussion of the nature of HD50138 

The variability in the line profiles, as cited by Pogodin dl997l l, 
seems to be caused by shell phases and/or outburst events. As 
described in the Sect. l3.ll from comparison of our spectroscopic 
data, we saw that there is a sort of temporal evolution, in partic- 
ular an increase of a blue emission component and a redshift of 
the absorption ones. Based on this, we believe that a new ejec- 
tion of material took place prior to our Narval observations, i.e., 
before March 2007 with increasing emission later on (our sec- 
ond FEROS observation). This new shell phase would be the 
responsible for the changes seen in our data. 

Another possibility is that the material was ejected in our 
direction, possibly by a hot spot on the stellar surface, form- 
ing a kind of "one-armed spiral", as suggested for classical Be 
stars. If this was the case, this would mainly explain the blue- 
shifted emission. On the other hand, even without having com- 
plete knowledge of the wind contribution for all lines, we sug- 
gested that for some of them, there is no contamination from the 
wind and they have a pure photospheric origin. Based on this, 
an extended and expanding atmosphere or a new ring or enve- 
lope with a lower outflow velocity that started to expand some 
time in early 2007 (or slightly earlier) could explain the changes 
in the radial velocities seen in the photospheric lines (see Fig. [9] 
and TableQ}. 

The results of Bjorkman et al. (1998), claiming the exis- 
tence of intrinsic polarization due to electron scattering, indicate 
a non-spherically symmetric structure with a gaseous disk seen 
almost edge-on, in association with an optically thin dust en- 
velope. Our results described in Sect. 13. 2. 2l confirm that such a 
scenario for the circumstellar dust is possible. 

Based on the stellar parameters derived in Sect. 13.2.31 we 
can indicate the position of HD 50138 on the HR diagram and 
compare it to evolutionary tracks. This is shown in Fig. [TBI 
where we plot the evolutionary tracks from Schaller et al. ( 1992 ) 
in the top panel for stars at solar metallicity. From this plot, 
HD 50138 is a (5.0 + 0.5) M star, since still on, or just evolv- 
ing off, the main sequence. The lower panel of Fig. [16] shows 
its position compared to pre-main sequence evolutionary tracks 
from Bernasconi & Maeder (1996). The ZAMS, as well as the 
birthline are also included. This plot indicates a stellar mass of 
roughly (6.0+1 .0) M , but for a star that is approaching the main 
sequence. 

From the position in the HR diagram alone, it is not possible 
to distinguish between a pre-main sequence and a post-main se- 
quence evolutionary phase, however we can discuss some points 
about the possible nature of HD 50138. 

- Pre-main sequence nature: 

The luminosity and the mass derived by our analysis agree 
with those expected for an intermediate mass pre-main sequence 
star. HD 50138 shows similar spectral and photometric varia- 
tions as seen in young objects. However, if it is a young star 
close to the main sequence, associated to its high temperature, 
its maternal cloud would already be dispersed and the remaining 
dust would be located in the accretion disk. This would be con- 
firmed by the forbidden lines seen in our spectra, which are not 
blueshifted, as seen in young deeply embedded stars (Waters & 
Waelkens 1998; Hamann |1994| ). In addition, based on its rather 
high temperature, HD 50138 would be a Herbig B[e] star. 

On the other hand, the absence of any nebulosity around 
HD 50138 and its situation far from any star-forming region 
(Pogodin 1997) is against a Herbig B[e] classification. The pos- 
sible shell ejection phases are also against a young nature. In 
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Fig. 16. Position of HD 50138 on the HR diagram compared to 
post-main sequence (top panel) and pre-main sequence (bottom 
panel) evolutionary tracks. In the bottom panel, the dotted line 
defines the zero-age main sequence, the thick solid line defines 
the birthline. 



addition, the results of Bjorkman et al. ( 1998 ) show that the dust 
is probably distributed in an optically thin envelope and not in 
an accretion disk. 

- A Be star close to (or just at the turn-off from) the main se- 
quence: 

The classification as a Be star that is still on, or just at the 
turn-off from, the main sequence cannot be ruled out at this 
moment. However, two points are against the classification of 
HD 50138 as a classical Be star: the forbidden lines and the IR 
excess due to dust. The formation of dust in an intermediate- 
mass star close to the main sequence seems to be a key problem 
for this scenario. 

The high-density circumstellar matter and the gaseous ion- 
ized disk-like structure speak for rather high mass loss. This 
mass loss must have been much higher than in the case of clas- 
sical Be stars, forming disk-like or spiral arm structures, if the 
material was expelled from some hot spots on the stellar surface 
and expanding outwards. Because of the much higher initially 
released masses in equatorial direction, the density in the outer 
regions will have ideal conditions for the production of the for- 
bidden emission lines. This mass loss might be less than in the 
case of B[e] stars. Then, in agreement with Jaschek et al. ( 1993 ), 
we might claim that HD 50138 is "a transition object between 
Be and B[e ] stars " . 
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- Binary system: 

For such a star to have a large amount of dust, as seen in its 
SED, another possible scenario would be linked to a binary com- 
ponent, where the dust could be formed in a long-lived circumbi- 
nary disk or by wind collisions. Baines et al. d20061 l suggest that 
the (possible) separation of components would be at least 0."5. 
However, considering the possible distance of HD 50138, this 
separation would be very large to have mass exchange (250 AU 
considering the smallest angular separation), which could ex- 
plain the spectral and photometric variability of this object. This 
would also hamper the possible classification of this star as FS 
CMa star (Miroshnichenko |2007l l, since this group of stars is ten- 
tatively classified as objects close to or still on the main sequence 
in binary systems with mass exchange. 

On the other hand, since the radius of the FEROS aperture is 
1", in principle, there might be a chance that we observed both 
components at the same time. However, we could not find any 
spectral feature to confirm the existence of a secondary compo- 
nent. This result also agrees with Corporon & Lagrange ( 1999), 
who searched for TTauri companions in a sample of HAeBe can- 
didates, based on several diagnostic lines, e.g. from Cai, Fei and 
Lii. They did not find any evidence of binarity for HD 50138. 
However, an HAeBe star as a companion could not be com- 
pletely discarded. Thus, at this moment we can neither confirm 
nor exclude a possible binary nature for HD 50138. 

- Other classes of B[e] stars: 

Since HD 50138 shows the B[e] phenomenon, we should 
also discuss the possibility that it belongs to one of the other 
known classes of objects that present it (Lamers et al. [1998). 
However, based on the absence of typical nebular and symbiotic 
lines in its spectrum, we can discard a compact planetary nebula 
and a symbiotic nature for this object. In addition, the low lu- 
minosity, in combination with no indication of nitrogen enrich- 
ment, exclude a supergiant scenario for this star. However, since 
almost 50% of the objects with the B[e] phenomenon are con- 
sidered unclassified, we cannot exclude that HD 50138 might be 
an object that represents a link between Be and B[e] stars. 

5. Conclusions 

HD 50138 is a very curious star that displays strong spectro- 
scopic variations. Based on analysis of new high-resolution data, 
we present a detailed description of these variations. Our analy- 
sis of the photometric data suggests that HD 50138 is a B6-7 III- 
V star, whose luminosity was tentatively obtained from a careful 
study of the influence of the possible circumstellar extinction 
sources and based on the new Hipparcos distance. A new value 
for the color excess, E(B — V) = 0.08 + 0.01 mag, was derived. 
In addition, we suggest that a new-shell phase or the formation 
of one-armed spiral could have taken place before 2007. 

Based on our results, a pre-main sequence star or a transition 
object between Be and B[e] stars, close to or just at the turn- 
off from the main sequence, or a binary scenario, can be neither 
confirmed nor discarded; however, an observational campaign, 
based on photometry, to derive a detailed light-curve and high- 
resolution spectroscopy associated to a detailed analysis in terms 
of the line profile appearances and variations have to be per- 
formed to confirm the possible shell phases and the slow-down 
of this star. A careful interferometric analysis, associated to a 
SED modeling considering different scenarios for the circum- 
stellar dust is in progress (Borges Fernandes et al., in prepara- 



tion) and will certainly provide better constraints for the circum- 
stellar geometry and the nature of this curious star. 
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Table 1. Lines identified in our FEROS spectra obtained in 1999 and 2007. The observed central wavelength {A t s ), radial velocity (y ra d), equivalent 
width of emission {W em ) and/or absorption (W abs ) components and identification of the lines seen in each date are provided. 
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2 Because of the complexity of the profile, the equivalent width of H/3 (and also of Hy and US) was derived considering the total one, i.e., the 
equivalent width of the absorption component minus that one from the emission components. 
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0.08 




Fen (m74) 6238.4 


6248.6 


48.0 




0.05 


0.06 


6248.7 


52.8 


0.03 


0.05 


0.06 


Fen (m74) 6247.6 


6301.1 


38.1 






0.55 


6301.1 


38.1 






0.58 


[Oi] (mlF) 6300.3 


6319.0 


38.0 






0.12 


6318.9 


33.2 






0.19 


Mgi (m23) 6318.2 






















Mgi (m23) 6318.8 






















Fen 6318.0 


6348.6 


70.9 




0.54 




6347.5 


18.9 


0.05 


0.32 




Sin (m2) 6347.1 


6364.6 


37.7 






0.20 


6364.6 


37.7 






0.21 


[Oi] (mlF) 6363.8 


6372.9 


70.6 




0.35 




6371.7 


14.1 


0.04 


0.22 




Sin (m2) 6371.4 


6385.6 


86.9 






0.11 


6385.1 


63.5 






0.18 


Fen 6383.8 






















Fen (m203) 6387.6 












6417.8 


42.1 




0.04 




Fen (m74) 6416.9 


6457.7 


60.4 




0.06 


0.06 


6457.1 


51.1 


0.10 


0.07 


0.10 


Fen (m74) 6456.4 


6564.8 


91.4 






61.35 


6564.1 


59.4 






57.56 


Ha 6562.8 


6680.1 


85.4 




0.34 




6678.7 


22.5 


0.03 


0.12 


0.03 


He i (m46) 6678.2 


7067.0 


76.4 




0.10 




7066.1 


38.2 


0.02 


0.06 


0.04 


He i (mlO) 7065.2 






















He i (m 10) 7065.7 


7496.5 








0.13 


7497.1 








0.15 


Uid 


7515.0 








0.31 


7514.2 








0.36 


Uid 


7713.0 


50.6 






0.26 


7712.2 


54.5 


0.11 


0.01 


0.06 


Fen (m73) 7711.7 


7773.5 


57.9 


0.20 


1.61 




7774.4 


92.6 


0.25 


1.77 


0.41 


O i (ml) 7772.0 






















Oi (ml) 7774.2 






















Oi (ml) 7775.4 


8325.3 








0.20 












Uid 


8346.9 


32.4 






0.31 












HiP23 8346.0 


8360.0 


35.9 






0.31 












HiP22 8359.0 


8375.5 


35.8 






0.61 


8374.8 


10.8 






0.86 


HiP21 8374.5 


8393.5 


39.3 






0.80 


8392.4 


0.0 






1.07 


Hi P20 8392.4 


8414.2 


32.1 






0.75 


8413.7 


14.3 






1.33 


HiP19 8413.3 


8440.0 








0.76 


8435.8 








0.85 


Uid 


8448.1 


46.2 






4.01 


8447.1 


10.7 






3.70 


Oi(m4) 8446.8 


8468.8 


53.2 






1.08 


8468.6 


46.1 






1.40 


HiP17 8467.3 


8501.4 


-38.8 






2.35 


8500.9 


0.0 






2.29 


H1PI6 8502.5 






















Can (m2) 8498.0 


8598.8 


14.0 






1.85 


8598.2 


-7.0 






1.86 


Hi P14 8598.4 


8630.7 


52.2 






1.14 


8630.6 


48.7 






1.73 


Ni(m8) 8629.2 


8664.9 


-3.5 






2.79 


8664.3 


-24.2 






2.94 


HiP13 8665.0 






















Can (m2) 8662.1 


8684.0 


20.7 






0.47 












Ni (ml) 8683.4 


8751.8 


44.6 






2.24 


8751.9 


48.0 






2.27 


Hi P12 8750.5 



